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Abstract. One-dimensional hard rod gases are explicitly constructed as the 
limits of discrete systems: exclusion processes involving particles of arbitrary 
length. Those continuum many-body systems in general do not exhibit the 
same hydrodynamic properties as the underlying discrete models. Considering as 
examples a hard rod gas with additional long-range interaction and the generalized 
asymmetric exclusion process for extended particles (£-ASEP), it is shown how 
a correspondence between continuous and discrete systems must be established 
instead. This opens up a new possibility to exactly predict the hydrodynamic 
behaviour of this continuum system under Eulerian scaling by solving its discrete 
counterpart with analytical or numerical tools. As an illustration, simulations of 
the totally asymmetric exclusion process (£-TASEP) are compared to analytical 
solutions of the model and applied to the corresponding hard rod gas. The case 
of short-range interaction is treated separately. 
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1. Introduction and outline 

Modeling stochastic many-body systems is fundamental to the investigation of driven 
diffusive systems [1-5]. So far, only few models have been solved analytically. With 
growing computer power, numerical tools like Monte Carlo (MC) simulations of 
discrete lattice models steadily gain in importance. However, for many physical 
processes, the restriction to discrete space is artificial, and the question arises, to 
which extent discretized models lead to a valid description of experimental settings 
taking place in continuous space. In the present paper, one- dimensional continuum 
models are investigated as to their explicit construction from discrete systems. The 
continuum gases treated in the following are one-dimensional hard rod gases exposed 
to an external driving field. Two types of hard rod gases are presented: The first 
consists of hard rods which do not show any interaction besides exclusion. This 
continuum gas is obtained as the limit of the discrete generalized asymmetric exclusion 
process for extended particles (^-ASEP) [6]. Basic properties of the discrete model 
and its hydrodynamic (HD) limit are summarized in section |21 before carrying out 
the construction of its continuum limit explicitly in section 13.11 Sections 13.213.31 are 
devoted to the investigation of the HD properties of the hard rod gas which turn out 
very differently from those of its discrete origin, the £-ASEP. The general proceedings 
of sections |21 and |21 are repeated for the construction and investigation of a second 
type of hard rod gas in section^ hard rods with an additional long-range interaction. 
In subsections 14. 11431 such a continuum gas with long-range interaction is realized as 
the limit of an exclusion process with headway-dependent hopping rates. This leads 
to remarkable results about the nature of discrete and continuum systems. Analysis 
shows that the hydrodynamic properties of this continuum gas are exactly the same 
as for the discrete ^-ASEP with constant stochastic rates. In section 1^31 the common 
macroscopic features of the £-ASEP and the long-range interacting continuum hard 
rod gas concerned are examined for their physical, microscopic causes. The usefulness 
of the correspondence established between those two systems is illustrated in section 
15. 21 where results from lattice MC simulations are applied to the continuum hard rod 
gas. 

2. A discrete model: the generalized asymmetric exclusion process 
(£-ASEP) 

The £-ASEP is a stochastic lattice model for a system of extended particles which 
interact via exclusion and which are exposed to a driving field. The microscopic dy- 
namics of the model define hopping processes of the particles taking place between 
neighbouring lattice sites with certain stochastic rates. The parameter £ fixes the 
length of the particles on lattice scale. For the case £ = 1, the £-ASEP reduces to the 
well-known asymmetric exclusion process (ASEP) [5]. FigureHillustrates the hopping 
dynamics for a lattice of L{— 35) sites, k — 1 . . . L, which contains N{— 4) particles, 
each covering £{= 5) adjacent lattice sites. The position of each particle is by con- 
vention taken as the last site that the particle occupies towards the right handside 
of the chain, and is marked with a cross. A particle may change its position from a 
site k to k + 1 with a rate p provided that site A: -I- 1 is empty. Vice versa, it may 
hop from site /c to fc — 1 with a stochastic rate q provided that site A: — ^ is not occu- 
pied by any particle. In the thermodynamic limit the number N of particles and the 
number L of lattice sites tend to infinity while the particle density remains constant. 
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On a coarse-grained scale the lattice constant a is a comparably small quantity. In 
a continuum limit a — s- 0, the discrete site label k may be replaced by a continuous 
coordinate x. The evolution of the local particle density p{x,t) is determined by a 
continuity equation of the form [6] 

dtp = -a{p-q)d^[vp] ~ a^^^d^^v + 0{a^) . (1) 
where v is a normalized average particle velocity [6,7] 

It must be stressed that this equation is an approximation for a system with discrete 
dynamics and that it does not describe the motion of hard rods with dynamics de- 
fined in continuous space. The £-ASEP was originally introduced as a lattice model 
for protein synthesis in 1968 [8, 9] and has been further investigated in this context 
recently [10-13]. In those works, the extended particles, hopping stochastically along 
a chain, represent the motion of ribosomes along the codons of a m-RNA template. 
For this purpose the description of the system by a discrete model is appropriate be- 
cause the codons form a real biological discrete lattice. Seeking to apply the £-ASEP 
to model a greater class of physical experiments which take place in continuous space 
(e.g. the study of colloidal suspensions), a continuum hard rod gas is constructed as 
a limit of the ^-ASEP. 



3. Continuum hard rod gas 

3.1. Construction of a hard rod gas: rescaling the particle length 

The derivation of the hydrodynamic equation ^ of the £-ASEP involves a continuum 
limit, in which the lattice constant a approaches zero. This equation in general does 
not describe the dynamics of a gas of hard rods in continuous space. The physical 
extension of the ^-ASEP particles la vanishes on a macroscopic scale. The particles 
appear like point particles from a coarse-grained point of view. 

In order to observe rods of non-zero extension in the continuum limit of space, their 
length on lattice scale must be rescaled during the limiting procedure as: 



A is held constant for a — > and indicates the true length of the rods in physical 
units. The continuum limit is performed for a fixed number of rods of length A on a 
chain of length A (L ^ oo sites). As the number t of sites which are occupied by one 
particle diverges for a ^ 0, the particle density p on lattice scale tends to zero, while 
the coverage density p'^ = Ip and the volume particle density p = ^ remain constant. 
This hard rod limit of the £-ASEP is illustrated in figure |3^b),A^B. 

3.2. Hydrodynamic limit 

The hydrodynamic limit of the continuum hard rod system is realized in a second step, 
considering a big system at constant density (A — > oo,^ oo, ^ = const.). During a 
second coarse-graining process, the length A of the rods becomes a comparably small 
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quantity (compare figure|3| (b),B^C). The hydrodynamic equation of the continuum 
hard rod gas may be deduced as the hard rod hmit of the HD equation Q of the 
discrete system (in figure|31 (a),B^(b),C). As the discrete particle density p becomes 
zero in the hard rod hmit, the coverage density p'^ wiU be considered instead. In 
terms of p'^, and employing the abbreviations B — p — q and S — equation Q is 
transformed into: 



dtp" = -aBdx 
= —aBdxp'^ + I 



{p^f 



a^Sdxx 



1-p- 



+ 0(0^) 



B 



A(l-p-) 



Sd.,r. 



+ 0(a3) 



(ifp^ ^ 1) 
(4) 

The scaling (■ = ^ and the limit of small a have been used to expand the terms. 
Furthermore I0J is only valid for the case of non-maximal density, which has been 
assumed in the first step in order to be able to perform the approximation without 
yielding singularities. The case of maximal density will be discussed in section IT^ 
Under Eulcrian scaling [dt = adx] , all terms in the current of orders higher than a are 
negligible. The evolution of the coverage density is determined by 



dtp"{x,t) + Bdxp"{x,t) = 0, 



p^^l. 



(5) 



The asymptotic approach of the mass current as a function of the coverage density 
from a convex to a linear function is demonstrated in figure |21 The linearity of the 
current in the hard rod limit implies that all rods are moving with constant velocity 



op 



(6) 



into the direction of the drive. In contrast to the discrete £-ASEP, interaction effects 
disappear completely. 



3.3. Solutions of the £-TASEP in its hydrodynamic and hard rod limit 

Equation jSJl excludes the case of particle distributions, where the density assumes its 
maximum at some point. In this section, explicit solutions of the Eulcrian current of 
the £-ASEP will be constructed. Their asymptotic behaviour in the continuum hard 
rod limit is discussed in particular for the case of partly maximal initial density pro- 
files. For the sake of simplicity of notation the i'-ASEP shall be limited to its totally 
asymmetric case, the so-called ^-TASEP, whose hopping rates equal p — I and q = Q. 

On Euler scale, only the lowest order terms of are taken into account for the 
^-TASEP current 

. p{l ~ ip) 

'^P^= l-i£-l)p 
which obeys the continuity equation 

dtp{x,t)+dxj{p)^0. (8) 

This partial differential equation (PDF) governs the evolution of the average particle 
fraction p per lattice site. In order to obtain the picture of macroscopic rods, it is 
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convenient to replace p by the volume particle density p. Furthermore every term ia 
is replaced by the constant length A. Equation © is thus transformed into: 



where 



dtp{x,t) + dxj{p) = 



1 — \p + ap 

In the following, solutions u{x,t) of this PDE of the general form 

m(1 — Xu) 



(9) 
(10) 



fin) 
u{x, 0) 



\ — Xu -\- au 

= cp{x). (11) 

are constructed for certain initial profiles (l){x). In the limit a ^ the solutions are 
solutions of the continuum hard rod gas. 



First, choose initial conditions of the form: 



Lo 

Rq 



X < 

x > 0. 



(12) 



The initial densities Lq and Rq are fractions of the maximal density: 
1 



Lq = 



Rq 



kLX 

1 

knX 



kL, k^ G 'Si . 



(13) 



The evolution of initial profiles which have a discontinuity are of special interest be- 
cause they determine the possibility of the formation of stable shocks in the system. 
One therefore demands fc^ ^ ku- Solutions of problem pi|l with the initial condition 
l|12|l are so-called weak solutions which are piecewise smooth and which obey (|ll|l only 
locally. 

Next, take a look at some general properties of f{u). Taking into account the physical 
condition u < j for any solution of problem (|ll|l . / proves to be a convex function 
{f'{u) < 0). The convexity of / implies that the unique piecewise smooth solution of 
(|nil,(|T2l is given by [14]: 



u{x, t) 




X < 
X > 



/(^o)-/(flo) 
fiLo]-fl'R„) 



Lo-Ro 
X < f{LQ)t 
f'{LQ)t <X< f'{RQ)t 

X > f{RQ)t 



Lq < Rq 



Lq > Rq 



(14) 



where 



(/') 



(15) 
(16) 
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is obtained from inverting the first derivative of / and only admitting types of solutions 
where f'{u) — y = j as 



The evolution of the system depends on the initial conditions: An initial discontinu- 
ity evolves as a stable shock for Lq < Rq, while for Lq > Rq, it dissolves in time 
according to a rarefaction wave solution. From a microscopic viewpoint this means 
that for the former case fluctuations which arise near the discontinuity are driven 
towards it according to the Rankine-Hugoniot jump condition for shock stability [14] 
Vieft > Vshock > Vright for the local velocitics around the discontinuity. In the latter 
case, fluctuations are driven away from the point of the (initial) discontinuity in the 
density profile and tend to soften it [5]. 

Limiting case a ^ 

So far, the solutions of the ^-TASEP current on Euler scale have been considered for 
finite a. In order to apply the solutions to a system of hard rods, the limit a — > is 
carried out in the following. For volume particle densities smaller than the maximum 
(p < ■^), one expects stability of any initial profile in accordance to equation 
However, it still remains an open question how an initial profile of partly maximal 
density evolves. The specifications = 1 and kji = 1 are therefore listed separately 
in the following. 

Case 1 (Lo < Rq) : 

• For LqjRq J {kL,kji ^ 1) the speed at any point of any profile is the same. 
Especially for (fT^ : 



Thus the discontinuity is moving with the velocity v shock = 1- 

• For Lq < i?o = J (fci > kfj = 1) the slope of the characteristics of the PDF (|ll|l 
starting at a; > for i = approaches 0, which corresponds to an infinitely fast 
propagation of any fluctuations in the initial high density regime towards the left. 
As the velocity at the right hand-side of the shock position, Vright, tends to minus 
infinity, the Rankine-Hugoniot jump condition is still satisfied. Thus even in the 
case of maximum density any initial profile is stable. The apparently infinitely 
fast spreading of the discontinuity has to be interpreted on lattice scale where 
any hard rod occupies infinitely many sites which have to be covered each time 
the fluctuation moves. 

Case 2 (Lp > Rq): 

• The case Lq,Rq ^ j{kL,kR 7^ 1) is similar to case 1, as the characteristics 
approach parallelism again, and the initial profile moves with constant velocity. 
The construction of h{y) is not necessary, because there is no space-time area 
which is not covered by the characteristics. 

• For Rq < Lq = J [kji > k^ — 1), the asymptotic behaviour of the 
rarefaction solution lima^o h{y) has to be considered for the space-time wedge 




(17) 



lim/'( — )^l(fce{A:i,fcfi}). 

a^O fcA 



(18) 
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lima-^o /'(j) = linia^o - ^ = -oo<2/<l = linia^o /'(fe^)- Expanding h{y,a) 
in orders of a yields: 

h{y, a) = ].- ^ + 0(ai) . (19) 

Thus h approaches a constant function in the limit a ^ 0: 

\imhiy,a) = \. (20) 

a-»0 A 

This means that the rarefaction wave solution approaches a solution containing 
a sharp jump in the limit, and that the discontinuity remains stable in this case 
as well (compare fig. EJ. 

Summing this up, the predicted result, i.e. a linear current in the hard rod limit, has 
been confirmed for the case of profiles of less than maximal density and its validity 
could also be enlarged to the case of piecewise maximal initial density profiles. 



4. Hard rods with long-range interaction 

Long-range interaction is introduced in the following as a new feature in the framework 
of the discrete £-ASEP. It is shown, how a generalized form of the ^-ASEP with 
headway-dependent hopping rates arises naturally from a mapping between exclusion 
processes and a certain class of zero range processes (ZRP) [15, 16]. A hard rod gas 
with long-range interaction is constructed as the continuum hard rod limit of those 
processes. 



4.I. The £-ASEP with headway- dependent rates and zero range processes 

Like the £-ASEP, a zero range process describes the stochastic hopping of particles on 
a lattice. However there are important differences in its microscopic definition: Firstly, 
the particles have zero interaction range, i.e. there is no restriction of the number of 
particles which may be located at the same lattice site. Secondly, the hopping events 
in general can take place between any two sites i and i' of the ZRP lattice and depend 
on the occupation numbers of those sites. In the following, a ZRP will be considered 
where the particles move along the lattice via next-neighbour hopping with rates p Un 
to the left and qun to the right, where ujn defines the occupation number dependence 
of the corresponding hopping rate and p—q determines the overall hopping asymmetry. 
For the case of constant rates = w, it has been described in [6] how to construct a 
one-to-one mapping between ZRP and ^-ASEP explicitly. The basic idea is to replace 
ZRP-sites by ^-ASEP particles and ZRP-particles by £-ASEP holes (compare figure 
[SJ. The densities pk of the extended particles and Ci of the ZRP particles are related 
as follows: 

Pk = ^ (21) 
where the discrete coordinates of the two systems transform like 

(22) 
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Application of transformations H21|l and (|22|l to ZRP with occupation number 
dependent rates u;„ yields a system of extended particles with exclusion interaction 
and headway-dependent hopping rates. 



4-2. The zero range process: stationary state and hydrodynamic equation 

The basic properties of the ZRP shall be shortly summarized in this section. The 
stationary state of a zero range process on a lattice of N sites j = 1 ■ ■ ■ N with rates 
w^Q) which depend only on the occupation r]{j) of the site of departure, is known to 
factorize into a product measure [15]: 

F*(77(l), ry(2), • • • , ,7(iV)) - firjil))firji2)) ■ ■ ■ f{7j{N)) (23) 



where 



with 



/(n)^ f"^ ^ , n = 0,1,2,. •■ (24) 



J2 (25) 
=1 

n 

n — ' Po = l (26) 



n=l 



Pr 

denote the single site probabilities for occupancy by n particles and are represented 
by power series in the fugacity z. 

Introducing a non-equilibrium analogue of the partition function 

oc 

Z = Y,Pnz" (27) 

n=0 

the ZRP density c can be deduced from the fugacity z as 

c = z-^lnZ (28) 
dz 

In the HD limit, the fugacity and the density become smooth functions of a continuous 
space coordinate y. The hydrodynamic equation for the density evolution of the ZRP, 
derived from the master equation of the stochastic process, takes the form [6]: 

dtciy, t) = aBdyz{y, t) -f a^Sdyyziy, t) + Oia^) . (29) 

Calculating (|28|l for a certain choice of rates LUn yields a hydrodynamic equation in 
terms of the ZRP density c{y, t). 



4-3. Construction of a continuum hard rod gas with long-range interaction 

As outlined above, a class of ZRP with occupation-number dependent rates through 
a one-to-one mapping induces a class of exclusion processes where the hopping rates 
of the involved particles are a function of their headway. This mapping is carried out 
explicitly in the appendix. Its implications for the investigation of continuum systems 
become evident in the present section where a hard rod gas with long-range interac- 
tion is yielded as the continuum hard rod analogue of a ZRP with occupation-number 
dependent rates. 
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Applying substitutions of densities and space coordinates arising from (PT|l and 
to equation H29|l yields the hydrodynamic equation for extended particles with exclu- 
sion interaction and headway dependent hopping rates of the general form (compare 
appendix, equation HA.2P '): 

dtp + aBd.,[pv{p)] + a^Sd^Mp) = (30) 

where the velocity term v(p{x)) = z{c{p{x))) is identical to the transformed expression 
of the ZRP fugacity. The hard rod limit of H30|l in terms of the volume particle density 
p = ^ under Eulerian scaling is calculated as (compare section HT^ : 

dtp= -Bd4pz{c{ap))] (31) 

where 

c{p) = c{ap) ^-{\-\). 

a p 

The form of the current-density relation for the hard rods is therefore completely 
determined by the fugacity- density relation 1)28(1 of the corresponding zero range 
process, i.e. by the choice of the form of the microscopic ZRP hopping rates w„. 

Current- density relation: determining the stochastic rates 

Having once established equation (|31|l . the rates of the ZRP which determine the 
basic relations 1(271 128(1 may be manipulated such that the corresponding hard rod 
gas shows a certain type of coarse-grained current-density relation. A hard rod gas 
exhibiting hydrodynamic phenomena similar to the ones of the ^-ASEP, for instance 
stable shocks, would be of much interest. Therefore one seeks to achieve a ^-ASEP type 
nonlinear current-density relation from equation 1(31(1 as a transformation from a zero 
range process with appropriately adapted stochastic hopping rates. For this purpose, 
the ZRP-rates ujn must induce an interaction taking place on macroscopic distances 
between any two rods. A convex current of ^-ASEP type can only be obtained for 
rates that increase with the distance of the particles which in turn is proportional to 
the number of free lattice sites n in front of a particle (i.e. the rates decrease with the 
particle density). 

The most simple choice of non-constant rates would be the one where w„ increases 
linearly in n or, if scaling is already taken care of, linearly in na. Calculation of the 
fugacity z{c{ap)) — ac{ap) = i — however, shows that the resulting current j = pz 
is also linear (and decreasing) in p. Though rates of higher order in n, for instance 
uJn = an^ , lead to a nonlinear current, such a current is concave and has a pole at 
p — Q. Figure El shows how the qualitative properties of the A-rod current depend on 
the choice of the ZRP rates as functions of n. 

The desired current, i.e. a current which resembles the ^-ASEP current also in the 
limit £ ^ oo, should increase with the density for low system densities and it should 
decrease with p for high densities. Therefore rates are necessary that are nearly con- 
stant for big macroscopic distances r and nearly linear in n for small macroscopic 
distances r = na. A possible choice is: 
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The derivation of c{z) for such rates H32|) is straightforward according to the relations 
given in section 14.21 The power series in z constituting the partition function is 
evaluated as 



Z 



oo ^ / n 



= (!-»)- 



\m—l 

a + l 



a + 1 



m — 1 



The fugacity z may be given explicitly as a function of c. From 

d a+l z 

c = z— ln(Z) = 

az a 1 — z 



follows: 



z{c) 



a+l 



a + 1 + ac 



In the limit a — s- and for constant volume particle density p the ZRP density c 
diverges like i. The fugacity in this limit approaches: 

ai(i-A) 
lim z(c(d)) = lim ° '^1,1 

1 - Xp 



1-{X-I)p- 

The hydrodynamic equation 13 1() applied to interacting hard rods under Eulerian 
scaling transforms into: 

p{l - Xp) 



dtp = -dxB 



(33) 



i-{x-i)p\ ■ 

Remarkably, one recovers exactly the same type of current-density relation as derived 
for the £-ASEP before taking the limit £ ^ 00 of macroscopic rods of rescaled length 
X = £a. 



5. Continuum gases and their lattice counterparts 

After a so far mainly formal derivation, a rather phenomenological approach to the 
comparison of discrete and continuum systems is taken in the following. Results from 
the preceding sections are reviewed focusing on their physical meaning as well as their 
practical significance for future calculations of similar systems. 

5.1. Types of interaction: l-ASEP and hard rod gas 

The coarse-grained behaviour of the continuum hard rod system with long-range in- 
teraction corresponds to the one of a lattice model of extended particles with exclusion 
interaction. The current-density relations in the continuum limit of both systems are 
formally equivalent. In the following, it is illustrated why those two systems macro- 
scopically evolve in the same fashion while they are governed by two different kinds 
of microscopic interactions. The procedure of coarse-graining during the HD limit is 
reviewed in greater depth. 
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Long-range interaction 

First of all, it shall be demonstrated that a long-range type of interaction is necessary 
for a hard-rod system with macroscopic interactions. The limit of vanishing lattice 
constant a is compared for the £-ASEP and the hard rod gas. 

• When performing the HD limit of the .^-ASEP, the real distance between two 
particles shrinks but the number of lattice sites between them remains constant. 
An interaction between two particles takes place, whenever a hopping attempt is 
rejected because neighbouring sites are occupied. Due to the finite distance, two 
particles which are separated only by holes in between, in general interact within 
a finite time. 

• Letting a approach zero for the hard rod gas, both the length of the rods and the 
inter-particle distances are rescaled like The real distance between two rods 
remains constant in this case, but the number of vacant lattice sites in between 
two rods approaches infinity. 

During the second coarse-graining, the real distance becomes comparably small 
again, but still contains infinitely many sites of the original lattice. Dynamics 
which are only defined on the microscopic lattice scale therefore can not possibly 
lead to any interaction. Two rods would need an infinite amount of time to meet 
and the sheer exclusion interaction does not play any role for the evolution of the 
system. 

If interaction effects are desired, one must equip the rods with some kind of long- 
range interaction which makes them feel the presence of each other even when on 
lattice scale at an infinite distance. The hopping rates H32I) have been scaled to 
a form such that they are functions of the macroscopic distance r instead of the 
diverging number of free lattice sites. 

Specific form of interaction: ZRP rates ujr = ^q^r 

Having understood the scales of interaction, its specific form shall be investigated. The 
long-range interaction of the hard rod gas is determined completely by the stochastic 
hopping rates of the corresponding zero range process and the discrete exclusion pro- 
cess. In the following it is explained why a certain coarse-grained hard rod current, 
in particular a ^-ASEP type current, must arise as a consequence of those rates. 



Consider a ZRP with a general form of hopping rates uj{r) which depend on the 
inter-particle distance r of the corresponding hard rod system like 

'^W-TT^- (34) 
Those rates result in a current-density relation of the form: 

J = A^- (35) 
1 — (A — x)p 

The choice of the parameter x implicitly fixes a length scale, determining the curvature 
of the function uj{r). The probability for a rod to leave its original position is 
determined by comparing its headway r to the length r + x. Choosing x = X the 
ratio = in average measures the local hole density of the lattice system. 
While a particle on the discrete lattice only hops if it encounters a hole in front due 
to exclusion, this extra hopping condition is artificially imposed on the hopping rates 
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of the rods. The resulting mass current for a; = A is symmetric and formally identical 
to the ASEP-current: 

J— = \j = p'{l - n . (36) 

Choosing x ^ A, equation (|34|l yields a non-symmetric current-density relation, 
especially for x = 1 the ^-ASEP relation. 



This means that the length scale fixed by x in the continuous system takes over 
the role of the lattice spacing in the discrete model. The rods are equally long as 
this unit in the case a; = A for the rod system and the case £ = 1 for the £-ASEP, 
resulting in a symmetric current. They measure A times and £ times respectively the 
length unit for a; = 1 or ^ > 1. In the case of the rods their real length A is essential, 
and the current is a function of the volume density. The length parameter £ of the 
£-ASEP particles instead indicates the number of covered sites, therefore the formally 
equivalent current is a function of the density per lattice site. 



5.2. Lattice simulations of continuum gases 

MC simulations are a standard numerical tool for predictions of the evolution of 
discrete models. In the previous section, the correspondence of a continuous to 
a discrete system (the ^-ASEP) has been established. This means that one may 
now employ such numerical approaches as exact simulations of a continuum system. 
Figures [3 and |H1 result from Monte Carlo simulations of the 3-TASEP as a chosen 
specification (^ = 3, i3 = 1, 5 = 0) of the ^-ASEP. They illustrate the temporal 
evolution of a hard rod gas on a periodic lattice with the current-density relation (|35|l 

.... P(l - 3p) 

starting with initial profiles of the form 

«M^0)^{'<-^ (37, 

A general form of this PDE has been solved analytically in section lXSl The numerical 
approach presented here involves MC simulations performed for an ensemble of 1000 
lattices of length L = 2000 each. Figure [S] depicts different stages of the MC evolu- 
tion of this discrete system where the initial profile H37|l is realized by a configuration 
with occupation of sites . . . 1499 by 500 3-mers and vacant sites 1500 . . . 1999. As 
periodic boundary conditions are applied, the initial profile at t = is a step with 
two sharp discontinuities at and 1499. For t > 0, hopping takes place according to 
£-TASEP dynamics (a particle hops to the right with rate 1 provided the lattice site 
is empty). One observes that the left discontinuity remains stable until it is reached 
by the particles coming from the break-up of the right edge of the step. At t = 480 
the equilibration process towards the stationary uniform density profile p = ^ has 
begun. Snapshots were taken at four times in the simulation, every 160 MC steps 
starting at t = after the lattice initialization, in order to show the softening out of 
the right edge of the initial step. In contrast to the TASEP {£ = 1), there is no linear 
decrease in the density profile but it takes the form which was already proposed by 
the rarefaction solution of the PDE (|ll|l in section 13.31 For a better comparison of 
simulation data and theoretical predictions, the rarefaction curve (|17|l is overplotted 
as a dashed line for times i > in figure |S1 In its region of validity, linking high- 
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and low-density regime, it matches the MC data very well. The remaining fluctua- 
tions are due to the discreteness of the MC system and the small ensemble number 
chosen for the averaging process. Figure [7| is a snapshot at t = 320 from the MC 
evolution of a similar but non-maximal initial step profile under 3-TASEP dynamics 
as a stable shock. The left discontinuity moves as a stable shock with constant veloc- 
ity towards the right. Except for stochastic fluctuations, the shock is absolutely sharp. 

The simulations performed match the expected results and confirm theory. The MC 
simulations, originally a tool to predict the behaviour of the discrete £-(T)ASEP, also 
describe a continuum system of long-range interacting hard rods. The existence of sta- 
ble shock solutions for this hard rod gas is demonstrated in the simulation pictures. 



6. Summary and Conclusions 



The construction of continuum hard rod gases has been achieved by appropriate scal- 
ing of and performing limits on discrete lattice models. For the example of the i!-ASEP 
as an exclusion process (model A) and a generalized process with additional long-range 
interaction (model B), two steps were performed: Firstly, in the hard rod limit the 
length of the extended particles in discrete space was rescaled to be infinite on lat- 
tice scale, but to stay finite on a macroscopic scale of continuous space. Secondly 
the hydrodynamic limit of those hard rods was gained by another coarse-graining 
procedure. The continuous models achieved in this way do not resemble the hydro- 
dynamic properties of the discrete models from where they arise. More remarkably 
instead, the hydrodynamic hard rod limit of model B yields exactly the same coarse- 
grained current-density relation as the hydrodynamic limit of the discrete model A. 
All basic physical properties of the ^-ASEP (hydrodynamic behaviour and phase tran- 
sitions [6,10,11,17]) are transferable to a continuum system of hard rods with a certain 
long-range interaction. The methodology behind this correspondence may be exploited 
for many other types of (lattice) models and inter-particle interactions. Given an em- 
pirical current-density relation of some real physical system, the established formalism 
allows for the construction of a lattice model which exactly reflects the hydrodynamic 
properties of the continuum system. This opens up new possibilities for the investi- 
gation of continuum systems with analytical and in particular numerical tools. 
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Appendix A. Generalized mapping between ZRP and exclusion processes 



Some key steps of the calculation of the HD equation of a generalized ^-ASEP with 
any type of interaction which can be mapped onto a ZRP as described in section 
are given. One starts with equation (|29|l . In order to rewrite its r.h.s in terms of p, 
one first carries out the outer derivatives of z{c{y,t)), yielding 



dtc{y, t) 



3z 

aB — dyc{y,t)~ 



— id 



.(A.l) 
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Secondly one introduces a function G{p) which is substituted for z{c) = z{c{p)), 
leading to 

dtc{y,t)= - aB{p'^)dpG{p)dyc{y,t) 

+ a'S [2p''dpG{p) + p'dppGip)] {dyc{y,t)f 

-a^Sp'dpG{p){dyyc{y,t)) 

and thirdly replaces the derivatives of c{y, t) by explicit expressions in p(a;, t) resulting 
in 

Px , a {px)'^ 



2 p5 



dtc{y,t)= aB{p')dpG{p) 

+ a^S [2p^dpG{p) + p%pG{p)] 

Hp.? 



{Px? 



a'Sp'dpG{p) 



+ 



Repeating all three steps for the l.h.s. 



where % is determined by 



dp 



dx 



dt 



— = aBz{y,t)+a'^Sdyz{y,t) - —Bdyz{y,t) 

one obtains the generalized hydrodynamic equation for particles with exclusion 
interaction and headway-dependent hopping rates: 

dtp + aB [G{p) + pdpG{p)] Px + a^S [dppG{p)ipx? + dpG{p)pxx] = 0. 

Thus: 

dtp + aBdx[pG{p)] + a'^SdxMp) = (A.2) 

and 

3{p) = aBpG{p) + a''SdxG{p). (A.3) 

The above formally introduced function G{p) has the physical interpretation of the 
average particle velocity for a normalized hopping asymmetry ij) — q) = 1, as in lowest 
order on the Eulcr scale the following relation holds true: 



<p) = f^=Gip). 



(A.4) 
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Figure 1. The £-ASEP: extended particles cover (. lattice sites each and move 
via next-neighbour hopping with stochastic hopping rates p and q. 




Figure 2. The ^-ASEP mass current as a function of the coverage density p'^: 
Asymptotically (^ — > oo), the current-density relation becomes linear. 
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(b) A-rods 



L lattice sites with N particles, each covering i sites, 
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(B) HD limit: p, p'^, I = const., a — > 0, i, -/V — > oo 



chain of length A with N rods of length A. 
volume densities p= ^, p" = 

(B) hard rod limit: p, p^, A, A = const., a — » 0, £ 

(C) HD limit: A 0, A, iV ^ oo, f = const 
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Figure 3. Limiting procedures for £-ASEP and A-rods (A — » B: continuum limit; B — » C: HD limit for A-rods). 
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Figure 4. The rarefaction solutions of the coverage density of the i- 

TASEP for different values of 1 {£ = 2, 10, 100) are plotted at time t = 1 over 
a certain x-range close to the zero-density regime {p'^{x, 1) = for x > 1). The 
initial coverage density was chosen as a step profile (p'^ = 1 if x < and p'^ = if 
X > 0). The steepness of the curves increases with the length £ of the particles. 
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Figure 5. Concept of the mapping between ZRP and f-ASEP: ZRP lattice 
sites are turned into particles, ZRP particles are replaced by holes; the stochastic 
hopping rates p and q are interchanged. 
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Figure 6. Qualitative dependence of the current-density relation j{p) of the 
continuum hard rod gas in its HD limit on the form of the corresponding 
microscopic ZRP hopping rates uj{n). 
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Figure 7. Evolved density profile after 320 MC steps of an initial step profile of 
non-maximal density for the 3-TASEP on a periodic lattice of 1500 sites. 
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Figure 8. Density evolution of the {£ = 3)-TASEP: The coverage density is plotted on the j/-axis against the lattice site numbers ■ ■ ■ 1999 
of a periodic lattice on the x-axis. Snapshots from the simulation data arc depicted for every 160 MC steps (t = 0,160,320,480). The 
dashed line indicates the theoretical rarefaction wave solution of the density profile. For better graphical illustration it is plotted over a 
bigger range than the one where it is a valid solution. In its definition range it perfectly matches the simulation data. 



